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We establish the existence of a static, classically stable 
string solution in a region of parameters of the generic two- 
Higgs Standard Model. In an appropriate limit of parameters, 
the solution reduces to the well-known soliton of the 0(3) 
non-linear sigma model. 

The effective theory of electroweak interactions may 
contain more than the single Higgs doublet of the mini- 
mal Standard Model. Theoretical arguments in favour of 
an extended Higgs sector Jl|] include supersymmetry and 
string theory, as well as the possibility of electroweak 
baryogenesis Q . Extra Higgs scalars give rise to a richer 
spectrum of non-perturbative excitations, such as mem- 
brane defects (|] and new unstable sphalerons Q . In this 
letter we will establish the existence of a new string exci- 
tation in the two-Higgs standard model (2HSM), and we 
will argue that it is stable in a region of parameter space 
extending into weak coupling. 

A systematic way to search for such non-topological 
excitations has been outlined by two of us in refs. [Spl: 
one considers appropriate limits of parameters so that 
field space acquires non-trivial topology, thus allowing 
for topologically-stable solitons. These solitons will gen- 
erally continue to exist when the parameters are relaxed 
by a 'small amount', so that the degrees of freedom that 
were frozen or decoupled in the limit affect only slightly 
the dynamics. A possible exception to this rule occurs 
when the limiting soliton has zero modes other than its 
center-of-mass position. Small corrections may lift in this 
case the degeneracy, and either fix or completely desta- 
bilize the excitation. This is exemplified by the Belavin- 
Polyakov soliton |Q, whose arbitrary scale can be fixed 
by embedding it into a gauged linear sigma model Q . 

The purpose of the present letter is to show that one 
can embed this same soliton as a stable static string de- 
fect in the 2HSM. The key observation is that in an ap- 
propriate limit of parameters the only relevant dynami- 
cal field is the relative SU(2) phase of the two higgses, 
partially-constrained to lie on a two-sphere S 2 . Our 
string defects are characterized by the fact that the map- 
ping of the transverse two-dimensional plane onto S 2 has 
non-vanishing winding number. In contrast to the pre- 
viously discussed vortex strings || the defects described 
here carry no net electroweak flux, have no symmetry 



restoration in their core, and do not exist in the one- 
Higgs standard model. 
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The Lagrangian of the 2HSM is C = 
lY^Y^ + \D^H 1 \ 2 +\D^H 2 \ 2 -V{H 1 ,H 2 ),whereW^ = 
%Wf - 3 V W« - ge abc W*W^ V = d^Y v - d v Y u . The 
physical Z° and photon fields are = W^cosOw — 



Y^sinOw and A^ — W^sinOy/ 



Y^cosOw and tanOw 



g' jg. Both Higgs doublets have hypercharge equal to 
one, their covariant derivatives are defined by D^Hi( 2 ) — 



dft + § gT a W^ + ^g'Y fJ _ J i?i(2) , while their potential reads 
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This is the most general potential (!]] subject to the con- 
dition that both CP invariance and a discrete Z 2 sym- 
metry (Hi — > —Hi) are only broken softly. The softly 
broken Z 2 symmetry is there to suppress unacceptably 
large flavor-changing neutral currents. 

Assuming all Xi are positive, the minimum of the po- 
tential is, up to gauge transformations, at < Hi >= 
(0,wi/%/2), < H 2 >= (0,w 2 /\/2). The perturbative 
spectrum consists of the electroweak gauge bosons with 
masses = g 2 [v\ + ^ 2 )/4 and mz — mw/cos9w , 



a charged Higgs scalar H + with mass 



v 2 )/2, and three neutral scalars, one CP-odd (A ) and 
two CP-even (h° and H°) , whose masses depend on the 
value of £. The neutral mass matrix simplifies consider- 
ably for sin£ cos£ = 0. For £ = 7r/2 the neutral masses 



are m 



Mn + «i)/2, m 2 



C — y/v and m z H0 



C + y/v, where ( = (Aiwf + A 2 v|) + (A 3 + jA 6 )(i;? + v 2 ) 
and n = (Ar + A 3 - ±A 6 )V + (A 2 + A 3 - + 

2v 2 v 2 ((A 3 + ±A 6 ) 2 + A 3 A 6 - (A 3 - iA 6 )(Aj + A 2 ) - AiA 2 ) . 

For £ = the masses of A , h° and H° are given by the 
same expressions with A5 and Xq interchanged. To sim- 
plify the discussion we take £ = 7r/2, g' = A3 = A4 = 0, 
and drop the decoupled U(l) gauge field V*. We will 
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comment on these assumptions at the end. 

Since we are interested in classically-stable static 
string solutions, we choose the Wq = gauge and work 
with static field configurations, so that Gauss' law is au- 
tomatically satisfied. Furthermore, for strings stretching 
in the X3 direction we take all fields to be independent 
of 23 and put W£ = 0. It is easy to verify that sta- 
ble minima of the ensuing two-dimensional energy func- 
tional correspond to stable infinite-string excitations in 
the original model. Unlike the Nielsen-Olesen-type vor- 
tex strings, the finite-energy solutions of interest to us 
will have |-Hi(2)| everywhere. This fact allows us 
to use the radial representation of the two doublets, and 
write them as #1(2) = ^1(2) ^1(2) (0 > !)■ ^1(2) are t w0 
positive functions, while ^1(2) are two smooth SU(2) val- 
ued functions on the plane. The space of smooth maps 
from the two-dimensional plane into SU(2) = S 3 is topo- 
logically trivial, and there is no topological obstruction in 
using the remaining freedom of smooth transverse-space- 
dependent gauge transformations to set either U% or U2 
equal to the identity matrix. We will choose U\ = 1 so 
that the most general Higgs configuration takes the form 



H 1 =F 1 



and H2 — F2 exp(— iOn • r) 



(2) 



with n • n = 1 and t the Pauli matrices. 



In the naive limit Ai, A2, A5 — * 00 the F\ 1 F% and 
freeze at their vacuum values Wi/v2> V2/V2 and £ = 
ir/2 respectively. This restricts the Higgs target space 
to a two-sphere S 2 , parametrized by the unit vector n. 
Sending A6 — > removes the potential of n. In order to 
decouple the gauge field we take furthermore the limit 
g — > 0. This is a priori dangerous because it renders the 
gauge fixing (Q) singular. To remedy the situation we 
must also take v± — > 00 while keeping the gauge-boson 
mass mw ~ gvi finite. We will see explicitly later on 
that a finite gauge-boson mass is indeed necessary for 
our stability analysis. 

In the above limit of parameters one is left with an 
effective 0(3) non-linear er-model C ~ u| (9 M n) 2 describ- 
ing the dynamics of the unit-vector field n(a;). The a- 
model has well-known topological solitons in two spatial 
dimensions Q , characterized by an integer winding num- 
ber N. This is the number of times n(ai) winds around 
S 2 as its argument covers the plane transverse to the 

defect. It is convenient to parametrize S 2 through a 

3 _ 



stereographic projection: n 1 



= 2n/(i + \n\ 2 ),n 



(1 - |^| 2 )/(l + |^| 2 ). The N = 1 string, extending along 
the X3-axis, and with the boundary condition n 3 — > 1 at 
00 is then given by 



= pe la /(z - zq) 



(3) 



where z = x\ + 1x2- The position zq , scale p and angle 
a are arbitrary parameters corresponding to soliton zero 



modes. The value n 3 = 1 at 00 will be imposed by the 
requirement of finiteness of energy once, as we do below, 
we let A 6 7^ 0. 

Following the same steps as in ref. J(|, we want now 
to relax slightly the above limits on the parameters, and 
study the fate of the solution (j^). A simple scaling ar- 
gument shows that the potential terms tend to shrink 
the soliton to zero size, while the gauge interactions tend 
to blow it up. Our task will be to show that in a re- 
gion of parameter space the soliton is stabilized at some 
fixed scale p. Our semiclassical analysis looks at first 
sight incompatible with the strong scalar-coupling limit 
considered above. To resolve this apparent contradic- 
tion, we must redefine the limiting theory in terms of 
classically- relevant parameters as in ||^J. To this end, 
we rescale the Higgs and gauge fields by my/ / \f2\~i, and 
the transverse space coordinates Xj (j — 1, 2) by 1/mw- 
The energy per unit string length takes the form 



£ = 



'W 



2Ai 



d 2 j 



F 2 2 (sm 2 e{d t n) 2 + (^0) : 



+ (d^) 2 + (d t F 2 ) 2 + -{Ff - i() 2 + f(R 
+ ^F 2 F 2 cos 2 6 + y (F^n 3 sin© - v^) 2 



r,2\2 
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We have defined g = gj\J2\\, A m = \ m j\\ for m = 
2, 5, 6, v\ = V2 cos (3/g, V2 = V2sin(3/g, tan/3 = v 2 /vi — 
V2/V1, and the current Jf = n a diB + sin 2 Qe abc n h din c + 
sin O cos Qdin a . We have kept the same notation for the 
rescaled coordinates and fields, as well as for the field 
strength W?j which is defined as previously but with the 
gauge coupling replaced by g. 

Since we will keep my/ different from zero in the lim- 
iting theory, we can use it to fix the length and mass 
scales by setting from now on mw = 1. We may, fur- 
thermore, treat Ai as the semi-classical expansion pa- 
rameter, which we will take sufficiently small so as to 
justify our semiclassical treatment. This leaves us with 
six classically- relevant parameters: v\, t>2, A2, A5, Xq and 
g. The limiting theory of interest can now be defined 
more precisely as follows: 

vi, X 2 — > 00; v 2 , A5 fixed; A 6 w 2 ,<7 — (3 — > 0. (5) 

This looks different from our original naive limit, but has 
the same dynamical effects. The first set of conditions 
freeze indeed F\, F 2 and O to their vacuum values, so 
that the only dynamically-accessible Higgs degree of free- 
dom is the field n(x), which parametrizes a two-sphere of 
non-zero radius. The second set of conditions decouples 
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the gauge field, and ensures that the dynamics of n is 
described by the usual 0(3) non-linear cr-model, without 
any additional potential term. The energy functional 



'(0) 



in this limit, admits the Belavin-Polyakov soliton as a 
(marginally-stable) solution. 

Relaxing slightly the limit (||) will introduce a poten- 
tial, unfreeze the "heavy Higgs modes" Fx — v\ = fx, 
F 2 — V2 = fi and 9 — 7r/2 = 6, and couple weakly the 
gauge field W° to the "light Higgs mode" n. These effects 
can be summarized by a classical (non-local) effective- 
energy functional M, 



X 2 v 2 



-fvl I d 2 y f k {x)G ki {x 



y)j?(y) 



+ 



(7) 



where j£ = e abc n b dkn c , and Gki is the two-dimensional 
massive Green function satisfying (<5jfcA 
$jk)G k i = S(x-y)Sji, and given by 



djdk 



Gki — 



1 



(27T) 



d 2 P - 
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' PkPi c ip-(x~y) 



(8) 



The first correction in ((?]) is a potential term. The second 
and the third come from the exchange of a heavy (radial) 
mode f 2 and of a vector boson, respectively. All three 
are small compared to £ ( ) , at least in the range of scales 
p 1 that will interest us. 

A series of comments on the above energy functional 
are here in order. First, £ e g contains an infinite se- 
ries of terms which come from integrating out classically 
the heavy Higgs modes and the vector bosons. Inspec- 
tion of (|4|) shows that all higher tree diagrams contain 
extra heavy propagators and/or extra powers of weak 
couplings. We can ensure that such corrections are in- 
deed subleading by insisting, for instance, that all three 
heavy Higgs masses be comparable. Secondly, had we 
set mw = 0, the gauge-boson exchange diagram would 
have been IR divergent for configurations with a simple 
power decay at infinity like (^). A consistent stability 
analysis would thus be impossible in this case. Finally 
we should point out that Skyrme-like and potential terms 
have been used to stabilize the Belavin-Polyakov soliton 
in the past j|. The difference is that in our case the 
non-local and non-renormalizable effective action (^) was 
derived from the classical field equations of the (renor- 
malizable) 2HSM. 

Following |6), we can look for a stable solution of 
(0) by minimizing the correction terms with respect to 



the collective coordinates of the zeroth-order solution 
(^). Since £ c g is independent of the center-of-mass po- 
sition and of the U(l) angle a, only the size p of the 
soliton is important. The current jf for the Belavin- 



(6) Polyakov soliton is jf — eijdji 
2p{xx — ix2)/{r 2 + p 2 ) and 3 



where 4> l 
~2p 2 /{r 2 - 



f i(j) 2 
P 2 )- 



A 



straightforward calculation then gives: 



5£ c s{p) 



ap 



dz 



z\K 2 (z) + K 2 {z)) 
(z 2 /p 2 ) + 1 



(9) 



where Kq,K\ are the modified Bessel functions, 5£ c r = 
£ c ff — £( ) , and we have defined the constants 



a = X 6 v(/2g 2 di , b = 8/3A 2 <Tv| , c = 2ng 2 v^/X 1 . 
The shape of 8£ c g(p), for different values of the param- 
eters a and b, was analyzed numerically. The results 
are given in Figure 1. At every point below the thick 
curve the energy has a local minimum, corresponding to 
a classically-stable soliton of size p = p(a, b). 




FIG. 1. The region of stability of winding string solitons. 



The soliton size stays constant along the thin straight 
lines, as shown in the figure in units of 1/mw- The as- 
sumptions that justified our perturbative treatment, and 
in particular the condition p ~ 1, are satisfied in a large 
part of the stability region. The validity of our analysis 
was confirmed independently by a numerical minimiza- 
tion of the energy functional of the 2HSM. We use, in 
an obvious notation, the most general axially symmetric 
ansatz 



Hx = 



hx + ih% 



H 2 = i[ ) , W 3 = W, 3 e r + Wjpe^ 



W+ = -^e^dWrx + iW r2 )e r + (W^i + iW^e^) (10) 
v2 

where W+ = (W 1 - iW 2 )/y/2. The ^-dependence is 
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shown explicitly and the ten real unknown functions de- 
pend only on r. 

Figure 2 shows the profile of the solution for g = 0.2, 
v 2 /vi = 0.3, Ai = 0.5, A 2 = 20.0, A 5 = 1.0 and 
Ae = 0.0001 (corresponding to the point A of Figure 
1) and with tension £ = 31.24^™^. We did not plot 
neither h(r) which is essentially equal to 6.7728 every- 
where, nor h 2 , W r i, W<p2, which to this accuracy 
are zero. As promised and contrary to what happens for 
Nielsen-Olesen strings, the magnitudes of the Higgs fields 
are everywhere non zero. Also, ~ 7r/2 at all points. 
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FIG. 2. The profile of a stable string. 

Let us comment briefly on the parameters which were 
not relaxed from zero in this discussion. From the anal- 
ysis in |(| we expect that turning on g' will favour larger 
soliton radii, without affecting our qualitative conclu- 
sions. The U(l) gauge interactions may even suffice by 
themselves to stabilize the soliton against shrinking. Sim- 
ilarly, turning on A3 should not affect significantly the 
discussion. Turning on A4 on the other hand, to render 
H + massive, presents a technical difficulty, because the 
corresponding potential term, evaluated for the zeroth- 
order soliton, diverges. The divergence is due to the slow 
approach of n 3 to its vacuum value at large distances. 
Although we don't expect this effect to be physically- 
significant, our minimization procedure would have to 
be modified in this case and one must check numerically 
whether H + can be made heavy enough to comply with 
experimental bounds. It is on the other hand a welcome 
fact of our analysis that A5 need not be large for stable 
strings to exist, and that Xq is favoured to be small. Re- 
alistic values of the CP violation parameter in the Higgs 
sector, which is proportional to the difference A5 — Ag [yj , 
may be consistent with string stability even if £ cannot be 
relaxed to very small values. Thus, there is no a priori 



indication that the experimental limits on the parame- 
ters of the 2HSM ]lj] are incompatible with the existence 
of our defects, but this must be decided ultimately by a 
numerical analysis of the stability region [|l0| . The possi- 
ble accelerator and/or cosmological implications of these 
solitons may of course depend sensitively on these results. 

Finally, having an S 2 target space in the effective the- 
ory suggests that one should also search for stable local- 
ized solitons classified by the Hopf index 7T3(5' 2 ) pJJ] . The 
analytic methods of this paper are not however applica- 
ble because the zeroth order sigma model lagrangian does 
not admit such Hopf solutions. 
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